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SUMMARY 

The unsteady flow field around an oscillating cascade of flat plates with zero stagger was 
investigated by using a time-marching Euler code. This case had an exact solution based on 
linear theory and served as a model problem for studying pressure wave propagation in the 
numerical solution. The importance of using proper unsteady boundary conditions, grid reso- 
lution, and time step size was demonstrated for a moderate reduced frequency. Results show 
that an approximate nonreflecting boundary condition based on linear theory does a good job of 
minimizing reflections from the inflow and outflow boundaries and allows the placement of the 
boundaries to be closer to the airfoils than when reflective boundary conditions are used. 
Stretching the boundary to dampen the unsteady waves is another way to minimize reflections. 
Grid clustering near the plates captures the unsteady flow field better than when uniform grids 
are used as long as the Courant-Friedrichs-Levy (CFL) number is less than 1 for a sufficient 
portion of the grid. Finally, a solution based on an optimization of grid, CFL number, and 
boundary conditions shows good agreement with linear theory. 


INTRODUCTION 

The accurate prediction of the unsteady flow field in turbomachinery remains a critical chal- 
lenge in predicting noise and aeroelastic response. Advanced designs of ducted propellers for 
large commercial aircraft have relatively large shroud diameters and swept composite blades. 

The blades are expected to be more flexible than the blades used in current turbofans. Blade- 
row interaction effects will be important for both aeroelasticity and aeroacoustics. Analytical 
and experimental predictions of the unsteady aerodynamics associated with these designs are 
needed to ensure safe operation and lower noise. 

Although three-dimensional methods are needed for a complete simulation of the unsteady 
flows associated with swept blades, two-dimensional analyses of linear cascades are used in prac- 
tice. Aeroelastic and aeroacoustic analyses commonly model a single blade row responding to a 
gust or to blade vibration. A number of numerical and analytical methods are available for 
addressing the two-dimensional, oscillating cascade problem. Time-linearized analysis has been 
used the most because it performs reasonably well near design conditions and is inexpensive to 
compute. Also, exact specification of the inflow and outflow boundary conditions is possible. 

An assumption made with all of these methods is that the unsteady flow field is a small pertur- 
bation of the mean flow. Work in this category includes Smith (ref. 1), Whitehead (ref. 2), 
Verdon (ref. 3), Hall and Crawley (ref. 4), and Fang and Atassi (ref. 5). Time-marching, non- 
linear schemes are also available that solve the full-potential (ref. 6), Euler (refs. 7 to 12), and 
Navier-Stokes (refs. 8 and 11) equations. These methods are needed for flows where small- 
perturbation approaches are no longer valid. However, they require larger amounts of computer 
time than linearized methods. In addition, proper treatment of the inflow and outflow boundary 
conditions, which is essential for accuracy, is difficult. 



This study fundamentally investigated the unsteady flow field around an oscillating cascade 
by using a time-marching Euler code. The code solves the nonlinear Euler equations by using a 
high-resolution, wave-split scheme. Only linear solutions are sought for a simple cascade of flat 
plates. Although using an Euler code is overkill for this problem, an exact solution that is 
known from linear theory can be used as a baseline for code validation. Furthermore, fundamen- 
tal studies of pressure wave propagation using a time-marching computational fluid dynamics 
(CFD) code are needed to better understand the limitations of numerical solutions. Choosing 
the proper grid resolution, time step size, and boundary conditions is important for resolving the 
unsteady flow field. Modeling the linear flow for oscillating flat plates is a first step toward vali- 
dating a numerical solution to be used for aeroelasticity and aeroacoustics and must be done cor- 
rectly before solutions for more complex flow fields can be considered accurate. This is 
particularly true for moderate blade passing frequencies associated with forced-response problems 
in aeroelasticity and with blade row interaction problems in aeroacoustics. 

In the present work the numerical predictions of the unsteady aerodynamics of an unstag- 
gered, oscillating cascade of flat plates were studied in detail. An unstaggered cascade was used 
so that the computational grid was perfectly orthogonal. An oscillation frequency was chosen 
that simulated a representative blade passing frequency in an advanced ducted propeller design. 
The importance of using the proper inflow and outflow boundary conditions was demonstrated. 
The effects of the boundary distance from the cascade, grid clustering, and time step size on the 
solution were investigated. 


NOMENCLATURE 


A Roe matrix 

c chord length, speed of sound, or unsteady characteristic variable 

ACP pressure difference coefficient, (p_ - p + )/(pjFi \h t |fc) 


e , energy 
F,G flux vectors 

Aj amplitude of oscillation for plunging motions based on chord length 

9m { } imaginary part of { } 
k semichord reduced frequency, o>c/( 2V[) 

M Mach number 

p static pressure 

Q dependent variable vector 

/2» { } real part of { } 
f.T time 

u,v velocities in x,y directions, respectively 

V total velocity 

x,y spatial coordinates 

curvilinear coordinate directions 
p fluid density 
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a 


interblade phase angle 
« oscillation frequency 

Subscripts: 

1,2 conditions at inlet or exit 

+,- upper or lower surfaces on plate 

LJR left or right of an interface 


GOVERNING EQUATIONS AND NUMERICAL ALGORITHM 

The unsteady, two-dimensional Euler equations written in conservative differential form are 
used and are transformed from a Cartesian to a time-dependent curvilinear reference frame. 

This transformation process and complete details of the numerical algorithm are presented in 
references 13 to 17. A brief overview is given here. 

The transformed equations can be written in vector form as 


dQ dF dG 

— + + 

dx 35 drj 


( 1 ) 


where 


Q = Jq 


(2a) 


F = J(Z,q * 5 X / + Z y g) 


(2b) 


G = J(y\ t q + r ] x f + q y g) (2c) 


and 


q = [p, p u, pv, e] T 
f = [p «, p« 2 + p, p«v, u(e + p)] T 
g = [pv, p«v, pv 2 + p, v(e + p)] T 


(3a) 

(3b) 

(3c) 


and 

q./m , and n are the metrics and J is the Jacobian of transformation, 
c x y • a y 
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These equations are discretized and solved by using a finite volume method, where cell 
centers are denoted as t and j : 


3g * V =0 

dx A£ Aq 


( 4 ) 


With A£ = Ail = 1 (by definition), this becomes 

f - 1 8 . F - »/<=) 


( 5 ) 


where 


6.0 = (). +1/2 - 0 .- V 2 ( 6 ) 


The components of the dependent vector Q represent average values over a particular cell. 
However, a method is needed to allow these fluxes to be accurately represented at cell faces. As 
discussed in reference 13, the method used in the present effort is based on the one-dimensional 
approximate Riemann solver of Roe (ref. 18) at cell interfaces for each coordinate direction. 

The method uses an approximate equation to represent a quasilinear form of a locally one- 
dimensional conservation law: 


dq 

dt 




= 0 


( 7 ) 


where A(q L , q R ) is a constant matrix that is representative of local cell interface conditions and 
is constructed by using so-called “Roe averaged” variables. Determining the eigensystem of A 
and knowing that the change in dependent variables across an interface is proportional to the 
right eigenvectors allow first-order flux formulas to be constructed. This approach of extracting 
flow field information from characteristically dictated directions is commonly called flux differ- 
ence splitting (FDS) and is applicable to multidimensional space, provided the assumption is 
made that all wave propagation occurs normal to a particular cell interface. In order to provide 
higher order spatial accuracy, a corrective flux is appended to this first-order flux. In addition, 
in order to control dispersive errors commonly encountered with higher order schemes, so-called 
limiters are used to limit components of the interface flux that result in total-variation- 
diminishing (TVD) schemes. All solutions presented herein were obtained by using the basic 
algorithm developed in reference 13, which is third-order accurate in space and second-order 
accurate in time. 


GRID 

For the simple flat-plate cascades considered in this study the flow equations were solved on 
Cartesian H-grids. Uniformly spaced grids were generated algebraically. Grids with cells 
clustered near the plates or grids coarsened near the inflow and outflow boundaries were 
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generated by using a two-dimensional version of the IGB code developed by Beach (ref. 19) for 
turbomachinery. This grid generator uses a hyperbolic tangent function to determine the spac- 
ing upstream and downstream of the cascade and ensures smooth variation in grid spacing. 

Once a grid is generated for a single passage, it is stacked to form a cascade with multiple 
passages. 

All applications are for plunging motions of plates with an interblade phase angle o of 90°. 
When the direct-store method presented in reference 12 is used, this requires four blade passages 
The Euler code is blocked so that the passage-centered H-grids and the corresponding dependent 
variables are stored on a solid-state storage device (SSD) on the Cray computer; the computer 
retrieves the information as needed. Within one block the lower computational boundary con- 
tains the upper surface of one blade in the cascade, and the upper computational boundary con- 
tains the lower surface corresponding to the adjacent blade. The grid from each passage is 
deformed so that it follows the plunging motion near the plates and remains fixed near the cen- 
ter of the passage. This is done by using weighting functions and is described in more detail in 
reference 12. The code time-marches the grid and the flow solution for harmonic oscillations of 
the plates. 


BOUNDARY CONDITIONS 
Solid-Surface and Periodic Boundary Conditions 

The solid-surface boundary conditions used in the present study were zero-pressure-gradient 
conditions. Janus (ref. 20) and others at Mississippi State University have found that zero- 
pressure-gradient boundary conditions are sufficient if the grid near the surface is adequately 
resolved. The boundaries that extend upstream and downstream of the plates are called periodic 
boundaries for plates with harmonic motion and a specified interblade phase angle. The Euler 
code stores phantom cells on either side of these boundaries and uses injection of the dependent 
variables to specify the boundary conditions for a given passage. 


Steady Inflow and Outflow Boundary Conditions 

The far-field, steady boundary conditions are based on characteristic variables and assume a 
locally one-dimensional flow at the boundary. A summary of their derivation that is based on 
reference 20 is given here. 

For the sake of deriving the boundary conditions the Euler equations are written in their 
nonconservative form: 

• V 

.y; dq dq ,dq n 

— + a— + b— = 0 
dx dl dq 


The matrices a and b are determined through an eigenvalue analysis. Multiplying by P ^ 1 and 
neglecting the derivatives in the q direction gives 
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( 9 ) 


-l dq 
' dc 


p: l p.A t p : 1 ^ - o 

5 C ? s 


where A £ is a diagonal matrix containing the eigenvalues X £ ; and P £ and P^ 1 are the left and 
right eigenvectors, respectively. The characteristic vector is defined as 


W t 


Fit 


( 10 ) 


where P £ is derived in reference 20 so that the elements of the characteristic vector become 


= 


w £ = 


m 


|V5| 



V 

1 

Q. 


1 c ')\ 


M 


p-* 


c, 

1 /J 


( 1 la.) 


(lib) 


w £ = 


J 

\p |V£| 

i/5|V£| 

. Pl c l 


(«,» + Zy V ) 


(11c) 


w l ~ 




- («,« + W 

Pl C l ' 


(lid) 


The corresponding eigenvalues are 


*{. k l = l' + + ly v 


(12a, b) 


X\ = + 5 y v) + c|V$| (12c) 

^ = a x u * S y v) - c|V5| 

The implementation of the characteristic variable boundary conditions requires knowing if the 
boundary is an inflow or an outflow. Generally, this is found by computing the sign of X to 
determine the directions of the characteristics. However, for cascades the direction of the flow 
and the characteristics at the boundaries are known. Let the subscript a denote approaching a 
boundary, b denote being on the boundary, and l denote leaving the boundary. For a subsonic 
inflow the characteristics approaching the boundary are set equal to the characteristics on the 
boundary by using equation (11). 


(12d) 
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(13a) 




P " 2 

Cl 


2 

Cl 


pM\ 


p.c 


li 


+ (*x“ + 5,V) 


mm 


p,c 


11 


{t x u + l y v) 


p - 




P 2 
C ') 


MM + (C x * - 5 y V) 

Pl c l 


MM 

pi c i 


'.(«,« + V) 


(13b) 

( 13c ) 

(13d) 


Solving equations (13a) to (13d) while assuming that the metrics at the boundary are equal to 
the metrics at a and l , gives the. following information for the flow variables at the inflow 
boundary: 

l . .... 

• • j : • • 

■ ‘ •’ 5 ;; Pb = \{P a + P,-* Pl C l[M“* - “l) + M V a " v /)]} : " r ' ’ ' ( l4a ) 

if- Z ’ • • 7 '• " ' 

(V4b) ‘ 


(14c) 

v» = v a + 5,^ (Hd) 

Pi c i 

s' . 

With the nonconservative flow variables known, the conservative dependent variables (eq. (3a)) 
can be easily computed at the boundary. 

The flow variables at the outflow boundary are found in a similar manner. The characteris- 
tics on the boundary are set equal to the characteristics leaving the boundary. In this case the 
incoming characteristic is the upstream-running pressure wave, and therefore the pressure is set 
equal to the pressure at l . For internal steady flows the pressure is usually specified by the 
user. The nonconservative flow variables become 




Pb - Pa 
P » = Po + ~ 


U b = U a + Zx 


Pa 'Pb 

Pl C l 


Pb = Pi 


(15a) 
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t>b = Pa + 


(15b) 


Pb - Pa 


, p a Pb 

- U a + lx 

(15c) 

P 2 C 2 


It 

+ 

At 

* 

i 

«■ 

(15d) 

P 2 C 2 



These characteristic variable boundary conditions are valid only for steady flows when applied to 
cascades because the outflow is assumed to have uniform static pressure. A number of flow 
solvers today mistakenly use similar boundary conditions for unsteady flows. For comparison 
purposes these boundary conditions will be called reflective boundary conditions when they are 
used later for unsteady flows. A better approximation for the correct unsteady boundary condi- 
tions is examined next. 


Unsteady Inflow and Outflow Boundary Conditions 

For unsteady flows in turbomachinery the exact specification of the unsteady boundary 
conditions for nonlinear problems is not feasible. Approximate “nonreflecting” boundary condi- 
tions are developed by assuming that linear theory can be applied. A very useful publication by 
Giles (ref. 21) presents the formulation of boundary conditions for internal flows. Much of the 
original theoretical work on this topic was done by Kreiss (ref. 22). Giles has transformed the 
theory into a form that can be easily implemented in general Euler solvers for turbomachinery. 
The basic implementation of this formulation determines the steady flow by using the boundary 
conditions derived herein and then solves the time-linearized form of the Euler equations at the 
boundary to determine the perturbation flow variables in terms of the characteristic variables. 
This implementation allows time variations of static pressure at the outflow and reduces reflec- 
tions from the boundaries. A brief summary that is based on reference 21 is given next. 

For the cases considered in this paper, x = £, y = tj, and t = T, so that the primitive form 
of the Euler equations can be written as 


dU 

dt 


. dU D dU 
+ A — + B — = 0 
etc dy 


(16) 


where 


U = (p, u, v, p) T 


(17a) 
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,5V. 



u 

P 

0 

(T 

A = 

0 

u 

0 

1 

P 


0 

0 

u 

0 


P 

W 

0 



V 

0 

p 

0\ 


0 

V 

0 

0 

B = 

0 

0 

V 

1 

p 


P 

0 

VP 

v > 


(17b) 


(17c) 


For small perturbations from the steady flow and with higher order terms neglected, a linear 
form of the Euler equations can be written as 


3U' 

dt 


+ A 


du' 

dx 


+ B 


dir 

dy 


= 0 


(18) 


where U' are the perturbation variables given as 


U' = (p, u, v,p) T 


(17a) 


and AJi are constant coefficient matrices based on the siteady flow variables. The nondimen- 
sional form of these matrices is 


b 1 0 O' 
0 u 0 1 
0 0 u 0 

f) 1 0 Uj 


(18a) 


B 


v 0 1 0 ^ 
0 v 0 0 
0 0 v 1 


(0 0 1 vj 


(18b) 


Following a Fourier analysis of these equations, a set of eigenvectors representing an entropy 
wave, a vorticity wave, and upstream- or downstream-running pressure waves can be deter- 
mined. The details of this procedure are given in reference 21. By using these eigenvectors and 
assuming locally one-dimensional flow at the boundary, the characteristic variables for unsteady 
flows are written in terms of the perturbation variables as 
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( 19 ) 


C 1 


-1 0 0 1' 


6p 

C 2 


0 0 10 


bu 

C 3 


0 10 1 


bv 

w 


>0 -1 0 !/ 


bp 


where 8p, 8a, 8v, and 8 p are the perturbations from the steady flow. 

For a subsonic inflow the amplitudes of the incoming unsteady characteristics (Cj, c 2 , c 3 ) 
are set to zero and the outgoing characteristic (cj is computed by using equation (19). For sub- 
sonic outflow c 4 = 0 and the remaining characteristics are computed. Once the characteristics 
are known, the perturbation variables are found by using an inverse transform: 


6p 

6m 

8v 

bp J 


-1 0 


0 

0 

0 


J. 

2 

J_ 

2 

0 

J. 

2 


n 

2 

1 

2 
0 
1 
2 


( c \ 

c i 


\ C 4J 


(20) 


The primitive flow variables can be found by using the following relationships: 

p = p 0 + 8p 


(21a) 


u = U 0 + bu (21b) 


v = v 0 + 8v (21c) 

P = Po + & P ( 21d ) 


where p^, u 0 , and p 0 are determined from the steady solution. Equation (3a) is then used to 
determine the conservative flow variables at the boundaries. 

RESULTS AND DISCUSSION 

A cascade of flat plates was used for studying the unsteady flow field caused by small- 
amplitude plunging motions. The purpose of the plunging motions is to propagate pressure 
waves upstream and downstream so that the accuracy of the CFD code in resolving the unsteady 
aerodynamics can be studied. Results are presented for one flow condition and various numeri- 
cal conditions, such as inflow and outflow boundary conditions, grids, and time steps (CFL 
numbers). Although the blade vibration problem is usually associated with flutter analysis, the 
forcing frequencies considered were high enough to simulate forced-response problems. The small 
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disturbances introduced by the plate motion allow for comparisons with small-perturbation 
theory. 


Model Problem 

The flat-plate cascade geometry used in the present study is shown in figure 1(a). The stag- 
ger angle was 0° and the gap-to-chord ratio was 1. Uniform flow with Mj = 0.50 was used for 
the steady flow. Small plunging oscillations with h x = 0.001 and a = 90° were used for all of 
the cases presented. A semichord reduced frequency it of 4 was chosen to simulate a realistic 
rotor-stator frequency in advanced ducted propeller designs. According to linear theory this flow 
condition propagates single-mode pressure waves both upstream and downstream in a direction 
not aligned with the Cartesian grids. These conditions define a model problem with a known 
solution that can be used to study wave dynamics by using CFD. 

The procedure used in the time-marching Euler code was to first determine a steady solu- 
tion with no blade motion. For flat plates with zero incidence this solution reduced to uniform 
flow and could be specified analytically. The unsteady solutions were started with each plate 
moving according to the specified interblade phase angle. Unless specified otherwise, the maxi- 
mum CFL number used in each of the following solutions was 1. The solution was run until the 
unsteady flow field reached periodicity. For the present study the unsteady pressures on the 
plate were F ourier transformed into their real and imaginary parts after each cycle of oscillation. 
The solution was stopped when the coefficients from the first harmonic did not vary signifi- 
cantly. This usually took 3 to 20 cycles of plate oscillation. The solutions presented herein were 
done on a CRAY YMP computer and took anywhere from 20 minutes to 5 hours of CPU time, 
depending on the grid and the time step. 


Visualization of Effects of Boundary Conditions 

A set of runs was done to visually demonstrate pressure wave propagation and the impor- 
tance of using correct inflow and outflow unsteady boundary conditions. Consider a uniformly 
spaced grid (dx = dy = 0.05 chord length) with the inflow and outflow boundaries located six 
chord lengths away from the leading and trailing edges, respectively. The grid resolution is too 
coarse for accurate computation of surface pressures but should be sufficient for visual demon- 
stration of pressure wave propagation to the far field. Figures 1(b) to (f) show “snapshots” of 
pressure contours just after the plates were impulsively moved from uniform flow. The amount 
of time between snapshots was approximately constant and is arbitrarily denoted by A T, which 
gives five instantaneous pressure plots over enough time for the downstream-running wave to 
reach the boundary. The left block used the “reflective” boundary conditions and the right 
block used the nonreflective boundary conditions. Because the reflective boundary conditions 
enforce constant static pressure at the outflow boundary, the downstream-running wave was 
reflected (figs. 1(e) and (f)). As the solution continued, the downstream domain was contami- 
nated with an upstream-running wave that was nonphysical for the desired solution. With 
proper grid resolution this wave entered the cascade and affected the unsteady pressures on the 
plates. 

Figure 2 shows the unsteady pressure distributions on the plate after the solution had run 
for 20 plate oscillation cycles. Four solutions are shown; two used the uniform grid with the 
reflective and nonreflective boundary conditions shown in figure 1 and the remaining two 
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solutions used the reflective boundary conditions on grids that coarsened as they extended to the 
inflow and outflow boundaries. Coarsening the' grid introduced numerical dissipation into the 
solution that dampened the waves as they propagated. (This type of grid has been commonly 
used for isolated airfoil analysis.) The grid had been coarsened so that the resolution in the pas- 
sage stayed the same (dx = dy = 0.05 chord length), but the aspect ratio of the cells gradually 
changed to 3:1 (dx = 0.15, dy = 0.05) at the inflow and outflow boundaries in one case and to 
10:1 (dx = 0.50, dy = 0.05) in the other case. There was a significant difference between the 
reflective and nonreflective cases that used a uniform grid, as would be expected from the con- 
tour plots in figure 1. Also, the solutions that used the grid-coarsening technique and reflective 
boundary conditions tended toward those for the nonreflective boundary condition case that used 
a uniform grid. It is expected that grid coarsening will only work if the grid spacing variation is 
smooth because abrupt changes in the grid may cause internal reflections. 


Effect of Boundary Distance From Cascade 

The effect of inflow and outflow boundary distance from the cascade was investigated by 
using both the reflective and nonreflective boundary conditions. The upstream and downstream 
boundaries were placed 0.25, 0.5, 1.0, 1.5, and 2.0 chord lengths away from the leading and trail- 
ing edges of the plates, respectively. The grid density was refined to dx = dy = 0.02 chord 
length. Each case was run through 20 oscillation cycles to give the unsteady pressure predictions 
shown in .figures 3 and 4. The results for the reflective boundary conditions (fig. 3) never 
reached a solution that was independent of the boundary distance. In contrast, the results for 
the nonreflective boundary conditions (fig. 4) were nearly identical when the boundaries were 
1.5 and 2.0 xhord. lengths away from the cascade. Thus, using nonreflective boundary conditions , 
can reduce the computational domain. 

Results from linear theory (ref. l) are also plotted in figures 3 and 4. The nonreflective 
boundary condition results that are independent of the boundary distance are in general but not 
exact agreement with the theory (short dashed lines in fig. 4). Because the theory is exact for 
flat plates, the differences between the solutions are due to numerical error. As demonstrated in 
the next section, the discrepancies are due to grid resolution. For now, the effects of boundary 
distance for different inflow and outflow boundary conditions determined in this study should 
still be valid. , . 

Effect of Boundary Conditions With Grid Coarsening 

The results shown in figure 2 have already demonstrated that the pressure waves by 
coarsening the grid should render a solution insensitive to either reflective or nonreflective , 
boundary conditions. In order to verify this, another set of solutions was considered that used a 
finer grid near the plates which smoothly coarsened to the inflow and outflow boundaries. The 
same boundary distances of 6 chord lengths upstream and downstream were used. The grid 
density was dx = dy = 0.02 chord length in the passage and coarsened to dx = 0.50 and 
dy = 0.02 chord length at the inflow and outflow boundaries. The unsteady pressure distribu- 
tions are plotted in figure 5 for solutions that used both the reflective and nonreflective bound- 
ary conditions, along with the exact solution from linear theory. The solutions were identical 
and differed from the theory by the same amount found in figure 4 for the nonreflective bound- 
ary conditions, Therefore, a solution has been found that is not influenced by reflections from 
the inflow and outflow boundaries. 
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Effect of Time Step (CFL Number) 


All of the previous solutions were run with a maximum CFL number of 1. For cases with 
varying grid spacing this means that the smallest cell did not exceed a CFL number of 1 and 
the larger cells had a CFL number less than 1. Grids with uniform spacing had a CFL number 
of 1 for every cell. Specifying a maximum CFL number of 1 preserved, at least formally, the 
second-order temporal accuracy of the numerical solution. However, many unsteady solutions 
that use a time-marching code run higher CFL numbers to reduce the run times. When this is 
done, an assumption is made that the higher CFL numbers occur only at the smaller cells and 
do not significantly affect the global solution. The following solutions show how wrong this 
assumption can be if care is not taken to ensure that the CFL number is small enough to resolve 
the unsteady flow field. 

The effects of CFL number were studied on two types of grid, including a uniformly spaced 
grid and a grid with points clustered near the plates. The nonreflecting boundary conditions 
were used in all of the remaining results. The uniformly spaced grid (dx = dy = 0.02 chord 
length) was identical to the one used in the study of boundary distance effects with the 
boundary located 0.5 chord length away from the cascade (figs. 3 and 4). Unsteady pressure 
predictions for three CFL numbers are shown in figure 6. The results for a CFL number of 1 
are the same as those shown in figure 4. The good agreement with linear theory is fortuitous 
because the predictions for a boundary distance of 0.5 chord length were found to be insufficient. 
Increasing the CFL number to 2 and 3 introduced temporal error, as expected. Using a CFL 
number of 3 completely destroyed the character of the unsteady waveforms. For this case two 
solutions have been plotted after 20 and 40 cycles of plate oscillation to verify that the solution 
has reached periodicity. 

A similar study is shown in figure 7 for a grid that was clustered in the y direction. The 
spacing near the plates was dx = 0.05 , dy = 0.02 chord length and smoothly varied to a spacing 
of dx = dy = 0.05 chord length at the center of the passage. Hence, the maximum CFL number 
occurred near the plates. The boundary distance was 0.5 chord length upstream and down- 
stream of the cascade. From the predictions shown in figure 7, clustering the grid near the 
surface helped preserve the character of the pressure waves but still quickly lost accuracy for 
maximum CFL numbers of 2 and 3. Note that for a maximum CFL number of 3, the minimum 
CFL number is 0.4 on this grid. This means that a portion of the grid should convect distur- 
bances correctly near the center of the passage. However, in the region of the grid near the 
plates, which was where the unsteady disturbances were being generated, the correct magnitude 
and phase were not computed. This is the reason why, although the nature of the disturbance is 
captured with a clustered grid, the accuracy is poor for higher CFL numbers. As demonstrated 
in the next section, it is possible for highly clustered grids to do a reasonably good job of 
predicting the unsteady flow field provided that an adequate portion of the grid has a CFL 
number less than 1. 


An Optimal Solution 

Finally, using the studies presented herein, an optimal solution is presented in figure 8. A 
highly clustered grid near the leading and trailing edges (dx = 0 . 003 , dy = 0.001 chord length) 
was used. The 121- by 41-grid (* by y) was stretched 6 chord lengths upstream and down- 
stream of the cascade by using the grid-coarsening technique. The maximum CFL number was 
10, which occurred near the leading and trailing edges of the plates. The CFL number quickly 
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diminished to less than 1 away from the plates, leaving a large percentage of the grid with an 
adequate time step to resolve the unsteady flow. Comparison of the unsteady pressure distribu- 
tions with linear theory in figure 8 shows good agreement for most of the chord length. 


CONCLUSIONS 

A fundamental investigation of the unsteady flow field around an oscillating cascade of flat 
plates with zero stagger was done by using a time-marching Euler code. The study demon- 
strated the importance of using proper unsteady boundary conditions, grid resolution, and time 
step size to accurately model the unsteady flow field. Imposing constant static pressure across 
the outflow caused reflections from the downstream boundary and was shown by plotting pres- 
sure contours at various times. These reflections contaminated the unsteady pressure distribu- 
tions on the plates. A study of the effect of boundary distance showed that the unsteady plate 
pressures never reached a solution that was independent of boundary distance. Using approxi- 
mate nonreflecting boundary conditions based on linear theory improved the solutions and 
allowed £he boundaries to be placed closer to the cascade (approximately 1.5 to 2 chord lengths 
away from the cascade). Coarsening the grid by stretching the boundaries away from the cas- 
cade introduced numerical damping to the solutions and was also an effective way to obtain cor- 
rect unsteady surface pressures on the plates, although it was only valid for unsteady flows that 
are due. to plate motion. These solutions were insensitive to reflecting or nonreflecting boundary 
conditions. Several solutions were presented that varied the time step size on both a clustered 
grid and a uniformly spaced grid. The results showed that the unsteady flow field became signi- 
ficantly inaccurate as the Courant-Friedrichs-Levy (CFL) number exceeded 1 on a uniformly 
spaced grid. The clustered grid results did a better job of capturing the flow but were still inac- 
curate because fewer grid cells had a CFL number less than 1. Good agreement with linear 
theory was found by using nonreflective boundary conditions on a highly clustered grid that 
stretched 6 chord lengths upstream and downstream from the cascade. For this case the maxi- 
mum CFL number was 10 and quickly diminished to less than 1 away from the plates. 
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(a) Cascade geometry for model problem. 
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Figure 1 . — Cascade geometry and Instantaneous static pressures for model problem. 
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Figure 1 . — Concluded. 
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Figure 6.— Effects of time step on unsteady pressure distributions when a uniform grid is used. 
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Figure 7. — Effects of time step on unsteady pressure distributions when a clustered grid is used. 
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